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Abstract 

The two dimensional lattice Ginzburg-Landau hamiltonian is simulated numerically for different 
values of the coherence length £ in units of the lattice spacing a, a parameter which controls amplitude 
fluctuations. The phase diagram on the plane T — £ is measured. Amplitude fluctuations change 
dramatically the nature of the phase transition: for values of £/a ~ 1 , instead of the smooth 
Kosterlitz-Thouless transition there is a first order transition with a discontinuity in the vortex 
density v and a sharper drop in the helicity modulus T. Both observables v and T are analyzed in 
detail at the crossover region between first and second order which occurs for intermediate values of 

e/o. 

pacs: 74.20.De, 74.60.-w, 64.60.Cn 



1 Introduction 

Several condensed matter systems undergo phase transitions which are triggered by vortex-like defects. 
For instance, this is the case of superfluid helium, superconductors and the melting transition. In par- 
ticular, the discovery of high-temperature superconductors has boosted the interest in the behavior of 
vortex lines in the mixed state and opened a new area which regards the physics of vortices as a new 
state of matter The Ginzburg-Landau (G-L) model involving a complex field ip = \tp\ exp[i9] captures 
the essence of all the above condensed matter systems. 

A common assumption when studying vortex-like topological excitations (singular phase 9 config- 
urations) is that amplitude fluctuations can be neglected and that the only relevant degree of freedom 
is the phase angle 9. The "phase only" approximation is equivalent to consider the limit in which the 
quartic coefficient of the G-L hamiltonian u — > oo which freezes the amplitude or equivalently to take the 
coherence length £ = 0. On the lattice this is the well known XY model. Thanks to the work of Berezin- 
skii Q and Kosterlitz and Thouless || in the early 70's we have a fair understanding of the XY model 
in two dimensions. It turns out that a topological phase transition, the so called Kosterlitz-Thouless 
(KT) transition, takes place driven by the unbinding of vortex-antivortex pairs at a temperature Tkt- A 
crucial ingredient of the KT transition is the logarithmic interaction between vortices. For that reason in 
ref. the authors concluded that in a charged superfluid, due to screening, vortices will always unbind 
at nonzero temperature. Latter on, it was realized that vortices in sufficiently thin superconducting films 
indeed interact logarithmically B-H and it was pointed out that the KT phase transition might also 
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apply to thin-film superconductors. The 2 dimensional flux line lattice (FLL) melting transition was also 
regarded as based on a KT-type theory |6) . More recently, it was suggested that some features of the KT 
theory might be present in under-doped superconducting cuprates 0. 

However, the nature of the phase transition triggered by vortices in 2D systems still remains under 
discussion. By means of Monte Carlo simulations of the XY model with a modified nearest neighbor 
interaction it was shown that, depending on the value of an additional parameter, continuous as well 
first-order transitions take place ||-[jl(|. The existence of both kinds of phase transitions is in accordance 
with the richer structure of the 2D Coulomb gas found by Minnhagen and Wallin 1 1 1 using self-consistent 
renormalization group equations and with the tendency towards first-order transition which develops 



in the case of a strong disorder coupling constant 12 . The 2D FLL also still remains controversial 



theoretically as well as experimentally. Experimental works and numerical simulations favor a KT-likc 



transition in some cases and a discontinuous one in others 13 . A recent extensive Monte Carlo simulation 



found a first order transition at a temperature close to the estimated one assuming a KT melting transition 

As Bormann and Beck |15| pointed out, at T ~ T c the amplitude fluctuations may affect the critical 
behavior and cannot be neglected. Therefore, instead of the XY model -with fixed amplitude fields- it is 
worth investigating the effects of vortices when amplitude fluctuations are taken into account. Depending 
on the chosen parameterization, amplitude fluctuations are controlled by the quartic coefficient u of the 
G-L Hamiltonian or by the coherence length £. Hence, in this work we simulate the Ginzburg-Landau 
model (with an amplitude and a phase degree of freedom per point) on a lattice of spacing a and obtain 
its two-dimensional phase diagram in the plane (T, £). We will sec that the nature of the phase transition 
of the G-L model diverges dramatically from the KT when the parameter u is chosen sufficiently small 
-which in fact is equivalent to take the coherence length £ ~ a- and that this is connected with the 
appearance of a discontinuous jump in the number of vortices. 

This article is organized as follows. In section II we define the lattice model, the observables to 
be computed as well as describe the used Monte Carlo procedure. In section III we report the main 
numerical results. Section IV is devoted to conclusions and final comments. 



2 Description of the model, the observables and the simulation 

We will consider square lattices of size L and spacing a and denote the lattice sites by x and the lattice links 
by (x, [i) with /i ~ 1, 2. A straightforward discretization of the continuum Landau-Ginzburg Hamiltonian 
produces the expression Jig ): 

im c = /3« 2 ElE -^)7« 2 +rm 2 +um% w 

X /J— 1 

where the superscript c denotes the ordinary parameterization in the continuum theory, (3 = jr^pp, Tn is 
the effective mass of the carriers and the coefficients r and u are analytic functions of the temperature, 

with u > for stability. Instead of working with this parameterization of the Hamiltonian -as we did in 

i 

ref. ]l7|j - it is better to introduce a dimensionless order parameter: ip x = (jp^ 2 ipx = 1 1 where \ipao\ is 
conventionally used because ip approaches this value infinitely deep in the interior of the superconductor. 
Hence the Hamiltonian can be written as 



PH C = = ^-a 2 - ^) 2 /a 2 + i(l - fe| 2 ) 2 ], (2) 

X fl — 1 
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where the lattice temperature Tl and the coherence length £ are given by 

Ti, 2muksT a 2 2m\r \ ' 

in what follows we will denote the lattice temperature simply by T, suppressing the subscript L. In the 
limit of £/a = (or u = oo) the radial degree of freedom is frozen and this model -sometimes said to 
describe soft spins with non fixed amplitude- becomes the XY model which is said to describe hard spins 
with fixed amplitude. A more interesting and less well studied limit is just the opposite i.e. £/a ~ 1 (by 
small values of the u parameter) which corresponds to large amplitude fluctuations. 
We have simulated the Hamiltonian (|2|) using a Metropolis Monte Carlo standard algorithm. The 
calculations were performed using periodic boundary conditions (PBC). In order to increase the speed of 
the simulation we have discretized the 0(2) global symmetry group to a Z{N) and compared the results 
with previous runs carried out with the full 0(2) group in relatively small lattices. For the case of Z(60) 
we found no appreciable differences. Lattices with L = 10,20,24,32,40 and (in some cases) 64 were 
used. For L = 10 and 20 we thermalized with, usually, 20,000-40,000 sweeps and averaged over another 
60,000-100,000 sweeps. For L = 40 and 64 larger runs were performed, typically 50,000 sweeps were 
discarded for equilibration and averaged over 200,000 sweeps. We also performed some more extensive 
runs near T c for small values of u. The errors for the measured observables are estimated in a standard 
way by dividing measures in bins large enough to regard them as uncorrelated samples. 
The following quantities were measured: 

i) The vortex density v. 

The standard procedure to calculate the vorticity on each plaquette is by considering the quantity 

m = ^-([9i - O2U + [0 2 - 0sW + [03 - 9iU + [04 - 0iW), (4) 
Ztt 

where [a^-n- stands for a modulo 27r: = a + 2im, with n an integer such that a + 2nn G (— 71", 7r], 

hence m — n\i + 7123 + ^34 + ^41 • If m 7^ 0, there exists a vortex which is assigned to the object dual 
to the given plaquette. Hence in the case d = 2, *m, the dual of to, is assigned to the center of the 
original plaquette p. The vortex "charge" *m can take three values: 0, ±1 (the value ±2 has a negligible 
probability), v defined as: 

v = *m x \, (5) 

X 

serves as a measure of the vortex density. 

ii) The energy density e =< H > / L 2 and the specific heat cy. 

Both observables were computed to measure the order of the phase transition, cy was computed 
using several procedures: simply as the energy variance per site i.e. cy = (< H 2 > — < H > 2 )/L 2 , as 
the mean of the variance over B bins per site and as the temperature derivative of the energy density 

Cy = iim^r. 

iii) The helicity modulus T fl8|| 

r measures the phase-stiffness. For a spin system with PBC the helicity modulus measures the cost 
in free energy of imposing a "twist" equal to LS in the phase between two opposite boundaries of the 
system. T is obtained in general as a second order derivative of the free energy with respect to 5 -which 
can be regarded as a uniform statistical vector potential- evaluated for S — > 0. In such a way one gets the 
following expression |]l9|, which generalizes the one introduced in ref. [Q, to an order parameter with 
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amplitude as well as phase variations: 



r = E ' I iMI 4>i I cos (^ - *i) > ~f < [ E ' I & II ^ I Bin(^ - e 3 )} 2 >}, (6) 

<ij> <ij> 

where the primes denote that the sums are carried out over links along one of the 2 directions (x or y). 

3 Numerical results 

First, in Fig. |l| we report the phase structure of that model in the (T, £) plane. The intersection point 
with the temperature axis £/a = corresponds to the XY standard model. Around £/a ~ 0.8 the phase 
transition line, separating the low temperature ordered phase from the high temperature disordered phase, 
changes from second order to first order. In order to illustrate this change in nature we present histogram 
results for two temperatures, one above and one bellow T — 0.8 for sizes L = 10, 20 and 40. The double 
peak structure of the energy histogram corresponding to the 2 coexisting phases, characteristic of a first- 
order transition, is showed in Fig. ||-a) for £/a = 0.85. Both peaks remain fixed as L increases and the 
width of each of them clearly scale as \J\-tjj) = j-j due to ordinary non-critical fluctuations, as can be 
checked from Fig. §-b). On the other hand, for £ < 0.75 the peaks are much lower and wider, they move 
towards to an intermediate value of the energy as L increases and their width do not scale as (Fig. 
l-c)). 

The central role played by vortex excitations in triggering and determining the nature of the phase 
transition can be seen in Fig. ||-(a) where v is plotted vs. T for different values of £ and L = 20. For 
£/a = 1 we observe a sharp jump in the vortex density v (triangles up). As long as we decrease £ the 
jump becomes more smooth and moves to higher values of T c until for £/a = 0.1 (+ symbols) we get 
something very close to the KT behavior (circles). The increase in the density of vortices when amplitude 
fluctuations are large, which is in agreement with the analytical computations of ref. Jl5| |, is due basically 
to the fact that amplitude fluctuations decrease the energy of vortices enhancing vortex production. The 
same happens for the XY model with modified interaction [ jTo| ]; in fact, the shape modification can be 
straightforwardly connected to a core energy variation. It is interesting to note that when £ ~ a the 
transition occurs almost from vortex (no bound pairs of vortex-antivortex) to a plasma of vortices, (see 
Fig. g-b)). When measuring v very close to T c from below we found that for £ = 1, n pa i r = v/2 becomes 
except for a very narrow interval -from T = 0.115 to T = Tc = 0.12- in contrast to what happens in the 
case of the K-T transition where bound pairs exist for a relatively wide temperature range -from T = 0.7 
to T = T c ~ 1-. It is interesting to compare the value of v for different values of £ in all the cases for 
the same very small reduced temperature t = (T c — T)/T = 0.01. We found (for different lattice sizes L) 
that ripair^/a — 1] is one order of magnitude smaller than the n pa i r [XY] = n pa i r [t;/a — 0] for the K-T 

transition; specifically: n p (^~ l (t — 0.01) ~ 0.0046 compared with n P ^~°(t = 0.01) ~ 0.05. 

The scarcity of bound pairs of vortex-antivortex below Tc for the large fluctuations regime (£ ~ a) can 
be explained in terms of the behavior of their free energy F pa i r — E pa i r — TS pa i r at T ~ T c from below. 

Roughly, E pa i r — 2£" c , where E c is the vortex core energy, and S pa i r ~ ^(jt)- The three quantities 
S pa ir,E c and T c all decrease as £ increases making difficult to disentangle the "energetic" contribution 
to v, proportional to exp[—E pa i r /T c ], from the "entropic" contribution ~ l/£ 2 . Fig. |4j shows a plot of v 
vs. £ always measured at the corresponding reduced temperature t = 0.01 and fittings with the energetic 
and entropic behaviors. For the intermediate range of £, although it is not easy to predict the behavior 
of the energetic factor, the entropy seems to be the main responsible for lowering the density of bound 
pairs. From £/a ~ 1, T c starts to decrease with A faster than E pa i r and thus exp[— F pa i r /T] decreases 
more and more sharply making smaller and smaller the probability of bound pairs of vortex-antivortex. 
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Figure 1: Phase diagram in the (T, |) plane for L = 40. The phase transition from the ordered phase to 
the disordered phase change from second order (dashed line) to first order (filled line) at £/a ~ 0.8 
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Figure 2: Histograms of e for L = 10 (below) ,L = 20 (middle) and L = 40 (above), a) £/a = 0.85: 
the 2 peaks becomes sharper and their position remain fixed as L increases, (b) Zoom of the right peak 
showing that its width scales as 1/L. (c) £/a = 0.75: the width of the 2 peaks do not scales as 1/L and 
they approach each other as L increases. 
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Figure 3: (a) v vs. T for £/a = 1 (triangles up), £/a = 0.8 (triangles down), £/a = 0.5 (x), £ = 0.1 (+), 
and the XY model(circles). (b) Zoom of ||-(a) showing the discontinuity of v at T c when £/a = 1. 



* : measured v (with error bars), o : ~1/s 2 , filled line: ~exp(-1/Tc} 
0.06 i 1 1 1 1 1 




Figure 4: v vs £/a (asteriscs with error bars) measured on a 20 x 20 lattice always at the reduced 
temperature t = 0.01 for that value of f/a. "entropic fitting" ~ l/^ 2 (circles) and "energetic fitting" 
~ exp[-E pair /T c ] (solid line). 
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Figure 5: e, v and r for £/a = 1 for sizes: L = 10 (circles), L = 20 (diamonds) and L = 40 (squares). 
Error bars for e and v are smaller than the symbol sizes. 
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Figure 6: e, v and T for ^/a = 0.5 for sizes: L = 10 (circles), L = 20 (diamonds) and L = 40 (squares). 
Error bars for e and w are smaller than the symbol sizes. 
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Figure 7: e, w and V for £/a = 0.1 for sizes: L = 10 (circles), L — 20 (diamonds) and L = 40 (squares). 
Error bars for e and v are smaller than the symbol sizes. 



Figures |||7] show plots of e, v and T for £/a = 1, 0.5 and 0.1. For £/a = 1 the transition is clearly 
first-order: we observe latent heat and discontinuous changes in v and r at T = T c . On the other hand, 
for £/a = 0.1 the results are similar to those of the KT transition. In particular, the energy is very close 
to the XY energy. 



4 Conclusions 

Therefore, in the G-L model the nature of the phase transition depends on the value of the coherence 
length parameter £ /a (which controls the thermal fluctuations of vortex cores) and the following picture 
emerges: 

1) For | > 0.8 the density of vortices experiments a discontinuous jump which coincides with a first 
order transition with large latent heat and a sharp jump in T all at a T c which decreases with £. This is 
the LGW-regime in which a sudden proliferation of vortices takes place instead of the unbinding of the 
much more smooth KT-regime. 

2) For | -C 1 we get basically the XY model (| ip \— 1) and the more subtle KT transition (with 
an unobservable essential singularity in the specific heat at T c and a much more small non-universal 
maximum above T c and a universal jump in T). The number of vortices and the energy evolve smoothly 
across the transition. 

3) For intermediate values of £ we have an interpolating regime between LGW and KT. 
Whether or not a large enough increment of * or u to alter the nature of the phase transition driven 

by vortices can be accomplished by varying some thermodynamic parameter, for instance the pressure, 
is something which deserves investigation. 

This change in the nature of the phase transition when amplitude fluctuations of ip are non negligible 
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(£/a ~ 1) is in agreement with very recent variational computations for the same model (g0[. Furthermore, 
we found first order transitions for £/a > 0.8, an interval which is included into the less restrictive of 
ref. Jl5| in which for £/a > 0.5 the authors found that the RG trajectories cross the first order line of 
Minnhagen's generic phase diagram for the two-dimensional Coulomb gas. 

Finally, our analysis could shed light on the nature of the melting transition. It is possible that 
the order of the melting transition in 2D depends on the particular conditions and details of the studied 
specimen which in turn translate into considerable different values of the coherence length 
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References 

[1] G. Blatter, M.V. Feigelman, V.B. Geshkenbein, A.I. Larkin and D. Vinokur, Rev. Mod. Phys. 66, 
1125 (1994) and references therein. 

V.L. Berezinskii, Sov. Phys. JETP 34, 610 (1972). 

J.M. Kosterlitz and D.J. Thoulcss, J. Phys. C6, 1181 (1973) 

M. R. Beasley, J.E. Mooij and T.P. Orlando, Phys. Rev. Lett. 42, 1165 (1979). 

A.F. Hebard and A.T. Fiory, Phys. Rev. Lett. 50, 1603 (1983). 

S. Doniach and BA. Huberman, Phys. Rev. Lett. 42 , 1169 (1979). 

J. Corson, R. Mallozzi, J. Orenstein, J.N. Eckstein and I. Bozovic, Nature 398, 221 (1999). 
A. Jonsson, P. Minnhagen and M. Nylen, Phys. Rev. Lett. 70, 1327 (1993). 
E. Domany, M. Schick and R. Swendsen, Phys. Rev. Lett. 52, 1535 (1984). 
J. E. van Himbergen, Phys. Rev. Lett. 53, 5 (1984). 
P. Minnagen and M. Wallin, Phys. Rev. B36, 5620 (1987). 
S. E. Korshunov, Phys. Rev. B46, 6615 (1992). 
K.J. Strandburg, Rev. Mod. Phys. 60, 161 (1988). 
Y. Kato and N. Nagaosa, Phys. Rev. B47, 2932 (1993). 
D. Bormann and H. Beck, Jour. Stat. Phys. 76, 361 (1994). 
M. N. Barber, Phys. Rep. 59, 375 (1980). 



G. Alvarez and H. Fort, cond-mat/0006341. In this letter we considered the usual parameterization 



in lattice quantum field theory for the Euclidean action. 

M.E. Fisher, M.N. Barber and D. Jasnow, Phys. Rev. A8, 1111 (1973). 

An analogous expression for the case of a granular superconductor appears in C. Ebner and D. Stroud, 
Phys. Rev. B28, 5053 (1983). In a similar way is calculated the orbital magnetic susceptibility, see 
for instance: A. Sewer, H. Beck, X. Zotos; Physica C 317-318 , 475 (1999). 

H. Beck and P. Curty private communication. 

Ch. Leeman et al, Phys. Rev. Lett. 64, 3082 (1990). 



9 



